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1 Introduction 



Finite temperature field theory is of great importance in many areas of physics. In the 
realm of high energy physics there are two main areas: phase transitions in the early 
universe and heavy ion scattering, where the phase structure of QCD may potentially be 
studied. One of the most salient features of finite temperature physics is that the effective 
degrees of freedom may be quite different at high temperature than at low temperature. 
Not the least of these changes is the potential change in effective dimensionality of a 
system at "high" temperature. We put "high" in quotes to emphasize that what is 
meant by high is a quite subtle issue, as are the nature and circumstances of dimensional 
reduction. Naturally, the question arises as to "high" relative to what? 

Many of the standard tools of field theory have been applied to finite temperature field 
theory, ranging from straightforward perturbation theory and other expansions, such as 
1/N, to various forms of the renormalization group (RG). The fact that understanding the 
high temperature regime is not trivial can be illustrated by some of the erroneous results 
that have been derived, such as that there is asymptotic freedom at high temperature in 
QCD, or that Ay? 4 at finite temperature exhibits a first order transition in the absence of 
an ordering field. 

As is well known, at temperatures T of the order of M /g or higher, the perturbative 
expansion breaks down. Here M = M(T = 0) is the zero-temperature mass and g is a 
generic zero-temperature coupling constant, for example g = ^/Ao (with Ao = A(T = 0)) 
in a scalar Ay? 4 theory. The reason for the failure of straightforward perturbation theory 
is the appearance of a thermal mass ~ gT which acts as an IR cutoff for the momentum 
integrals in Feynman diagrams and generates inverse powers of g in the perturbative 
expansion. The physical origin of these infrared problems is that at high temperature the 
effective degrees of freedom are quite different to those at zero temperature. 

Through a careful analysis of perturbation theory it can be shown that these prob- 
lems can be cured by an appropriate resummation of the most infrared divergent graphs 
at high temperatures, the so-called hard thermal loops (HTL) [1]. In the case of Ay? 4 
theory, it is sufficient to replace the zero-temperature mass M in the propagators by 
the temperature-dependent mass M(T), calculated consistently order by order, while the 
coupling constants may be treated perturbatively once the mass corrections have been 
resummed (see [2] for the explicit results to two loops in the resummed theory). 

However, the HTL approximation is invalid near a second order or weakly first order 
transition where the temperature-dependent mass M(T) becomes small and the resulting 
IR divergencies are of a qualitatively different nature: the infrared fluctuations dominate 
in the same sense as in critical phenomena, where they invalidate a perturbative treatment 
of all physical quantities. For instance, the existence of another RG fixed point invalidates 
a perturbative expansion around the Gaussian fixed point which in this case is infrared 
unstable. A key property of the physics near the phase transition is universality - 
universal quantities near the phase transition are completely independent of the zero 
temperature parameters of the theory. This is clearly not the case in the regime which can 
be described by the HTL resummation, where the zero temperature coupling is affected 
only slightly by thermal fluctuations. 

Henceforth, we will refer to the regime near a phase transition, with M(T) <C M , gT, 
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as the universal high temperature (UHT) regime, in contradistinction to the non-universal 
high temperature (NHT) regime where T ~ M /g or larger but M(T) is not small in the 
above sense. A typical example of a relevant UHT situation would be electroweak theory 
in the vicinity of the electroweak phase transition (assuming that the phase transition 
really is, at most weakly of first order), while an example of a NHT situation would be 
the electric sector of QCD at high temperatures. One of the main theoretical tools to 
access the NHT regime in QCD has been precisely the HTL expansion [1]. The question 
then is: under what physical circumstances is one regime reached versus the other as a 
function of temperature and to what extent the two regimes may be treated on an equal 
footing in a given calculational scheme? 

The answer to the question of which regime is reached by increasing temperature 
depends on whether or not a state of broken symmetry exists at T = 0. If the state at 
T = is symmetric then the mass M(T) tends to monotonically increase as a function 
of temperature. At temperatures gT ~ M , the perturbative expansion around the zero 
temperature theory breaks down in a non-universal way (NHT regime). Eventually, a 
regime is reached where gT ^> M and the physics becomes independent of the value of 
the zero-temperature mass. 

If, on the other hand, there is a broken symmetry at T = 0, characteristically the mass 
will diminish until a phase transition is reached. If the phase transition is second order 
or weakly first order then near the phase transition M(T) <C Mq and perturbation theory 
will again break down, but now in a universal way (UHT regime). However, there may 
also exist a NHT regime below the transition, where gT ~ M but M(T) is not (yet) very 
small relative to M , hence in this regime perturbation theory will break down in a similar 
manner to the case without symmetry breaking. Above the transition temperature, in 
the symmetric phase, the mass will tend to increase thus entering (again) a NHT regime. 
For M(T) ^> Mq, the physics is expected to be insensitive even to the symmetry of the 
zero temperature state. 

Realizing that the regime of small mass M{T) is dominated by universal fluctuations 
attempts have been made to describe it using techniques gleaned from critical phenomena. 
An example of such an approach can be found in [3] . An alternative approach has been to 
describe the regime numerically [4] . Such attempts have all been based on an assumption 
- that a dimensional reduction takes place and that the low energy limit of the theory 
may be described by an effective three dimensional theory. In the case of [4] the parameters 
of this effective theory are calculated in terms of the parameters of the four dimensional 
theory and the temperature. However, the matching of the four dimensional to the three 
dimensional theory is done via a perturbative analysis analogous to the HTL expansion, 
so the obvious question arises: is a dimensionally reduced NHT regime reached, in which 
case a perturbative matching can be performed, before having to fully take into account 
all important universal infrared corrections? (Note that in the case of a non-abelian gauge 
theory a perturbative treatment is only possible below the critical temperature, i.e., in 
the phase of broken symmetry.) These questions in fact, and answers to them, are the 
principle subject of this paper. 

As far as we are aware the only techniques capable of unambiguously answering 
these questions have their origin in the RG and, in particular, must be "environmentally 
friendly" [5] RGs in that they must be capable of describing quantitatively the change 
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in effective degrees of freedom in the system as a function of "environment" - - in this 
case temperature. It naturally follows that such RGs must necessarily be temperature 
dependent. However, temperature dependence is not sufficient to guarantee success as 
various temperature dependent RGs have been used without success (see for example [6]) 
to access the UHT regime. The important requirement is that the renormalized finite 
temperature parameters give an adequate description of the physics at any temperature. 

Such environmentally friendly RGs have been used to successfully investigate both the 
NHT and UHT regimes within one calculational scheme, both running the thermal mass 
and the temperature. For instance, in [7] the critical temperature and amplitude ratios 
were calculated for a Ay? 4 theory, thus showing that the critical temperature, albeit a non- 
universal quantity, could be calculated directly using RG techniques. The emphasis in 
this paper is on using environmentally friendly renormalization to understand when and 
where a dimensionally reduced description of the system is appropriate, what is its nature 
- universal or non-universal, and where and when important calculational techniques 
such as the HTL expansion are valid. We take, again, a scalar theory as a testing ground. 

The structure of the paper is as follows. In Section 2 we briefly present our finite 
temperature renormalization prescription and the resultant flow equations. In Section 3 
we present the solutions of the flow equations in the NHT and UHT regimes, and consider 
the validity of dimensional reduction and the HTL expansion and finally in section 4 we 
draw some conclusions. 



2 Finite— temperature renormalization 

We define our model via the bare Euclidean action at temperature T — 1/(5 



S[<p B ] = J^dtJ d d ~ l x 



\iy VB f + l -MW B + ^ B 



(i) 



In d — 4 ultraviolet divergences require renormalization of the parameters \ B and M B , and 
also a wave-function renormalization. In d = 3 the theory is super-renormalizable with 
respect to ultraviolet divergences. However, the massless theory can give rise to infra- 
red divergences. Hence, we will define renormalized parameters via the normalization 
conditions 



^T ( t 2 \p,^T),M(r),\(r),T = r) 



p=0 



= 1, 



T^( P = 0^(t),M(t),X(t),T = t) = M 2 (r), 
rf ) (p = 0,v?(r),M(r),A(r),T = r) = A(r). 



(2) 
(3) 
(4) 



The renormalization is at an arbitrary temperature r, hence the renormalized parame- 
ters are temperature dependent and therefore, in principle, "environmentally friendly". 
Note that the normalization conditions (2) and (4) are applied to the "transverse" vertex 
functions and which are defined via the conditions 
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which is just the equation of state for the model in the presence of an external current J, 
<p is the renormalized field expectation value, and 



, \ r {4) 

p(2) = p(2) + U_ -2_ 



(6) 



Although the transverse functions are motivated by considerations of the O(N) model, 
where (6) is a Ward identity for instance, they are equally valid in the analytically contin- 
ued limit N — > 1 and lead to significant advantages in the renormalization when compared 
to normalization conditions on the "longitudinal" vertex funtions. More details can be 
found in [7]. Note that the physical mass m(T) at the normalization point is given by 



m 2 (r) 



M 2 (r) 



1 + r{ 4) (r) 



(7) 



p=0 



and coincides with M(r) when tp — but not otherwise. 
The RG flow equations are 



T ~dr ~ 7v ' 



T- 



dM 2 (r) 
dr 



a. 



Mi 



_d\(r) 
dr 



(8) 



and take different functional forms in the broken and symmetric phases. In a perturbative 
treatment of these flow equations we treat A in perturbation theory, but for M in each 
diagram we eliminate M\ in favour of M 2 with the aid of condition (3). This eliminates 
all diagrams which contain a tadpole as a sub-diagram. Subsequently, we differentiate 
with respect to r and solve for the flow functions, expanding them to the loop order we 
are working. It has been verified to two-loops that the resulting flow equations are free 
of UV and IR divergences. 

The beta functions at one loop are 



l v = 0, 
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M 



r a d®_ 

2 T dr ' 
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d 



r>T r 



, r <T C , 



* = -2 A V 



(9) 
(10) 

(11) 



where the symbol © stands for the one-loop diagram with k propagators, without vertex 
factors, at zero external momentum. It can be obtained from the following basic diagram 
in d dimensions (in the present case, we will consider the limit d — > 4) 



M 



O = r £ fJL-JL M e + (2 m ,f + 

n=—co \ J 

T(-f)M d 2r d r«> 

~ (47r) (d - 1)/2 r(^) Jo 7? 



(Air) 
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where z = M/t, by differentiations with respect to M 2 . The first derivative gives 



<l = ®, (13) 



dM 2 



whereas for k > 1 we have the general rule that the derivative with respect to M 2 of 
the loop with k propagators gives (—k) times the loop with k + 1 propagators. The 
integration of equations (10-11) was used in [7] to investigate the critical regime, in 
particular showing how the critical temperature could be derived using RG methods and 
also calculating critical amplitudes associated with the transition. 
Equation (11) is easily solved, giving 

\-\r) = A-Vo) + \ [<2> (M(r), r) - @ (M(r ), r )] , (14) 

where to ensure that the same initial conditions are imposed on both sides of T c one may 
use the requirement that the bare coupling is the same in both phases to find 

\-\t + ) = AlV-) + \ [<2> (M + , r + ) - (2) (M_, r_)] . (15) 

After solving the flow equations and setting the arbitrary RG temperature r equal to the 
physical temperature T the parameters M(T) and A(T) describe the behaviour of the 
vertex functions and at zero momentum. 

Here, we have used temperature as an RG flow parameter. There is also benefit in 
using a fiducial value of the finite temperature mass as flow parameter. In this case 
defining the floating coupling h = 4AM 2 (3) , chosen so that the quadratic term of the beta 
function has unit coefficient, one finds 

8h 

M-^ = ~(A-d cS )h + h 2 (16) 

where d e g defines the effective dimension of the system [5]. By changing variable M — > 
z = M/t we have d e g(z). In the limit z — > oo one finds that d e s — > 4, with h being 
proportional to the zero-temperature coupling, while in the limit z — ■> one finds that 
d c ff — > 3. Thus, we see that d e s can usefully be employed as a measure of the effective 
dimension of the system. Note that this definition of effective dimension is universal. 
In other words universal quantities, such as critical exponents, will interpolate between 
those associated with different dimensionalities. d e s will not a priori be suitable as a 
measure of non-universal dimensional reduction as in that case all "universal" quantities 
are mean-field like due to the suppression of infrared fluctuations. 



3 Results and Comparison with other Methods 

In this section we will present results obtained by integration of the differential equations 
(10-11) for the temperature dependence of the mass and coupling constant. We considered 
both the case when the zero temperature theory exhibited a broken symmetry and when 
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it did not. Additionally, in both these cases we considered different initial conditions for 
the zero temperature coupling. 

As mentioned, our principle concerns are to examine the phenomenon of dimensional 
reduction and also the regime of validity of the HTL approximation. Normally dimen- 
sional reduction is something that is assumed based on simple dimensional arguments. 
For example, Fourier transforming in the imaginary time formalism one finds an infinite 
tower of masses M5 2 + 47r 2 T 2 n 2 . One argues that in the infrared only the lowest mode, 
n = 0, is relevant, the higher mass modes decoupling in the infrared for sufficiently high 
temperature T (T ^> Mo), thus leading to an effective three dimensional theory. The 
problem with this is, as we have already discussed, there are two quite distinct regimes 
where this can take place — the UHT and NHT regimes — and, also as we have em- 
phasized, these are universal and non-universal respectively. Hence, in the UHT regime 
dimensional reduction will be a universal phenomenon, independent of the zero tempera- 
ture theory, whereas in the NHT limit its existence will depend on the zero temperature 
parameters. In the UHT regime, due to the existence of universal fluctuations, one would 
expect to see three dimensional critical behaviour with associated three-dimensional crit- 
ical exponents. This is indeed what is found. However, the NHT regime is not fluctuation 
dominated due to the suppression of infrared fluctuations by the large thermal mass. In 
this case one expects the behaviour to be "mean field" like. Once again, this is what is 
observed. 

On a qualitative level, in the imaginary-time formalism, /3 — 1/T is the extension of 
space-time in the temporal direction. As long as the correlation length £ = 1/M is much 
smaller than (3, or equivalently T <C M(T), the system does not "feel" the restriction 
in the temporal direction, while for ^ > or T > M(T), the field configuration can 
be considered as approximately constant in this direction as long as we consider modes 
with four-momenta p < T, and the system reduces to an effectively three-dimensional 
one. This criterion replaces the naive one according to which dimensional reduction 
occurs whenever T 3> M from a consideration of the bare propagator in the Matsubara 
sum representation. The difference consists precisely in considering the renormalized or 
"resummed" mass M(T) instead of the zero-temperature mass M . Thus, we may fix 
a criterion for dimensional reduction that there is dimensional reduction when z(T) = 
M{T)/T < z , where z is some suitably chosen number. Typically, we will choose 
z = 0.3 so that dimensional reduction is considered to be valid when the system "size" 
is 0.3 times the correlation length. In the plots of M(T) vs. T below, a fixed value of z$ 
corresponds to a straight line through the origin with slope Zq. For the part of the curve 
M(T) below this line we then consider a dimensional reduction to have occurred while 
above it not. 

Of course, the crossover from four to three dimensions is smooth. The effective di- 
mension introduced in the previous section takes this crossover into account in a natural 
fashion. For instance, using the effective dimension, if we consider the regime where 
z(T) = 0.3 we may ask to what extent there is a dimensional reduction, i.e. to what 
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Figure 1: The effective dimension defined in Eq. (17), as a function of z — M/T. 
extent the system appears to be three dimensional. Explicitly, 

d cS = 4 + ^" M ) 2 = 6-6M>® . (17) 
dhiM 

which turns out only to depend on z. We note that it is essential to keep the logarithmic 
term in the expressions for the one-loop diagrams because it represents the zero temper- 
ature contribution. We have plotted d e s as a function of z in Fig. 1. The limit z — > oo 
corresponds to r — > (for a finite mass Mo), where physics is entirely four-dimensional, 
whereas z — > holds at the phase transition where M(r) — > at finite r, and a universal 
dimensional reduction to three dimensions is fully realized. We can use d e & to justify our 
choice of z Q = 0.3 noting that for this value of z the effective dimension is equal to three 
with good precision (less than 3.01). Further evidence for what values of z$ lead to dimen- 
sional reduction can be gleaned from considerations of any universal physical quantity. 
In particular in [5, 8] effective critical exponents that interpolate between two different 
dimensions were considered in the context of a finite size system in a film geometry of 
thickness L with periodic boundary conditions (mathematically equivalent to the present 
case). It was found that the crossover region is broader crossing over from four to three 
dimensions rather than three to two. Also, it was found that the crossover took place at 
smaller values of z$ for four to three crossover than for three to two. Taken together these 
two points led to the observation that three-dimensional critical behavior was not realized 
until T ~ 100M(T) at least. In light of this evidence our criterion of z = 0.3 is quite 
liberal as, although in this case there were deviations from four-dimensional behavior, 
one could not argue that the behavior was effectively three-dimensional. 

Our other principal concern is the validity of the HTL expansion. The HTL expansion 
devises a systematic scheme of partial resummation of the standard perturbation series 
in the case of finite temperature, thus recovering an expansion in powers of the coupling 
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constant (actually, in powers of y/Xo in the present case), while a naive perturbative 
expansion runs into problems for "soft" external momenta p < y/\)T for temperatures 
of the order M /\Ao or higher, because the thermal mass of order V^o^, acting as an 
infrared cutoff, generates inverse powers of the coupling constant in the loop integrals. 
Parwani [2] showed the consistency of the HTL resummation up to two-loop order for 
the case of A0 4 theory. In the following, we will refer to his results generally as "the HTL 
approximation" . 

Up to one loop in the present scalar theory, the HTL resummation amounts merely to 
a resummation of daisy diagrams, which systematically incorporates the thermal mass to 
order V^Ao T into the bare propagator. However, in the UHT limit the HTL approximation 
does not properly take the thermal corrections to the coupling constant into account. 
These thermal fluctuations are crucial for a correct description of the phase transition. 
More generally, we will compare our expression for the coupling constant, 

A(r) = J , (18) 



^ + i <2/(*o-)) 



with the HTL one-loop result [2] 



Ahtl(t) = A -^©'(z(t)). (19) 

In Eq. (18), we have neglected the logarithmic term which becomes notable only at ex- 
tremely high temperatures, consequently (2) is replaced by (2) , defined as the negative 
of the second M 2 -derivative of O where the first, explicitly r-independent term is left 
out in Eq. (12). In the original one-loop HTL expansion, the mass appearing in the loop 
integral in Eq. (19) is replaced by its expansion to first order in V^Aoj a difference which 
we neglect here for the sake of easier comparison. We also include the contribution of 
the zero-temperature mass M in z(t) = M(r)/r, which is usually put to zero in the 
HTL approach. We consider the discrepancy between the two formulas (18) and (19) as 
appreciable as soon as 

Ahtl(t) 

where we typically will take v = 0.1, corresponding to a 10% difference between the two 
values, which is equivalent to 
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0.3. (21) 



The latter relation defines a value of z for any given Ao, which is the value represented by 
the HTL curves in the figures below. 

Our criterion (20) for the validity of the HTL approximation might be criticized as 
subjective. To improve on this matter, observe that there is precisely one value of the 
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Figure 2: The HTL-critical coupling constant \ v as a function of the parameter v speci- 
fying the criterion according to Eq. (22). 

initial coupling constant Ao such that the equation (20) is exactly fulfilled in the limit of 
large r (still neglecting any logarithmic contribution). We call this value the HTL-critical 
coupling. Generalizing (20) to 



we can define an HTL-critical coupling \ v for every value of v. For all initial coupling 
constants Ao < X u , the HTL approximation will then be acceptable according to the v- 
criterion, for sufficiently high temperatures. In Fig. 2, \ u is plotted as a function of u, for 



Having motivated our criteria for dimensional reduction and the validity of the HTL 
expansion we will now consider the explicit results: In all the figures, the temperatures 
where dimensional reduction (approximately) applies are the ones where the curve M(T) 
lies below the dimensional reduction curve (short dashed line), whereas the HTL resum- 
mation is adequate wherever M(T) lies above the HTL curve (dotted line). In Figs. 3-5 to 
the left, we see graphs of M(T) and A(T) as a function of T for the three different initial 
couplings 0.05, 0.5 and 2.0 respectively where we start at T = in the symmetric phase 
with M = 1. For small coupling, A = 0.05, we see that the HTL approximation is always 
valid and that there is a non-universal dimensional reduction starting at about T = 3Mq. 
Note that straightforward perturbation theory is expected to be applicable as long as 



VAoT <C M , hence in the present case for T <C 4M . For larger couplings A = 2 NHT 

dimensional reduction does not set in until about T = 5M and the HTL approximation 
is only valid for T < 6M . One immediate conclusion of this is that the existence and 
the value of the "high" temperature above which HTL is invalid is sensitively dependent 
on the initial coupling. Hence, contrary to current folklore, even for large temperatures, 
dimensional reduction and HTL only apply for sufficiently small initial values Ao, from 
Figs. 4, 5 roughly A < 2.0 and A < 1.0, respectively. The breakdown of the HTL ap- 
proximation is due to the fact that thermal corrections to the coupling are significant, as 




(22) 



< v < 1. 
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M =1 .0, X =0.05 M =1 .0, X =0.05 




5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 

T T 

Figure 3: Mass M (solid line) and coupling constant A (dashed line) as a function of 
the temperature T, as obtained from integrating differential equations (10-11) with the 
initial conditions M = M(T = 0) = 1.0 (in arbitrary units) and Ao = A(T = 0) = 0.05. 
For the plot to the left, the system is in the symmetric phase at T = 0, while for the 
plot to the right, it is in the phase of broken symmetry. The lines of constant z(T) = 
M(T)/T correspond to the criteria for dimensional reduction (short dashed line) and HTL 
approximation (dotted line) discussed in the text. 



M =1 .0, A. =0.5 M =1 .0, X =0.5 




Figure 4: As in Fig. 3, but with initial condition A = 0.5. 
M =1 .0, A. =2.0 M =1 .0, A, =2.0 




Figure 5: As in Fig. 3, but with initial condition A = 2.0. 
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at high temperatures they will vary as Aoa/Ao which when resummed gives rise to a large 
contribution. 

Note that if we are more rigorous about our criterion for validity of the HTL expansion, 
for example demanding agreement with the full thermal coupling to 5% rather than 10%, 
then the HTL line will be steeper and the critical coupling for which a full fluctuation 
analysis is required in the high-temperature limit correspondingly smaller. Also, if we 
require a stricter criterion for dimensional reduction then the corresponding dimensional 
reduction line would be flatter. Clearly, for a certain value of the zero-temperature 
coupling constant depending on the criteria chosen, the HTL and dimensional reduction 
lines cross thus indicating that the regime of dimensional reduction can no longer be 
described within the HTL approximation if the coupling constant is larger than this 
critical value. 

In Figs. 3-5 to the right, we see analogous results to those in the same figures to the left, 
but now for the case where we start in the broken symmetry phase with M = 1. In this 
case there exists a phase transition at a critical temperature T c . At T c , the mass as well as 
the effective coupling constant (when defined as in (4)) vanish. For large temperatures, the 
mass grows proportionally to the temperature and the coupling constant tends towards a 
constant value. For both quantities, the values at large temperatures become independent 
of the T = initial condition Mo, while they do vary with Ao- Notice in particular that the 
asymptotic value of A can deviate appreciably from Ao, at least if Ao is sufficiently large. In 
actual fact, at temperatures very much larger than the ones shown in Figs. 3-5, M(T)/T 
and A(T) show a logarithmic rise which originates in the usual renormalization of the UV 
divergence in the one-loop contribution to T^, and finally leads into a Landau pole. The 
high-temperature behaviour of the solutions will be analyzed in more detail below. In 
Figure 3, for A = 0.05, we see that the HTL approximation is valid except in a small 
region around the critical temperature, of width about 0.1T C , and that there is a non- 
universal dimensional reduction starting at about T = 0.2T C with a universal dimensional 
reduction setting in much closer to the critical temperature where universal quantities 
are essentially three-dimensional. On the other hand, in Figure 5, for A = 2.0, we see 
that non-universal dimensional reduction does not set in until T > 0.7T C and moreover 
that the HTL approximation breaks down at about the same temperature. Once again, 
the critical value of the coupling constant where the HTL and the dimensional reduction 
lines cross and hence the regime of dimensional reduction cannot be described within the 
HTL approximation, depends on the criteria chosen. The important conclusion to be 
drawn here is that there may be no non-universal dimensional reduction via which an 
effective three-dimensional theory may be posited before universal fluctuations invalidate 
the use of an HTL approximation necessary to analytically calculate the parameters of 
the three dimensional theory as is done in the context of the electroweak theory in [4]. 
Once again we emphasize that our present estimates are based on a liberal criterion for 
granting dimensional reduction and a fairly liberal one for determining the validity of the 
HTL approximation. 

Several properties of the solutions of the flow equations (10-11), and in particular 
the limit of high temperatures, become more transparent when one considers the flow 
in the z-A-plane. To begin with, we will neglect the logarithmic contributions which 
become important only at extremely high temperatures, thus considerably simplifying the 
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following discussion. The integration of the flow equation for A defines, after neglecting 
the logarithmic term, a characteristic curve X(z) along which a fixed physical system 
evolves with changing temperature. By considering the limit r — > or z — > oo we infer 
that the characteristic curve is specified by only choosing the initial coupling constant 
Ao- Furthermore, the matching condition (15) implies that the curve is the same below 
and above the critical temperature, for a fixed physical system (which is in the broken- 
symmetry phase at r = 0). 

We now determine the flow of the system along its characteristic curve for increasing 
temperature. The flow equations for the mass imply that 

dz A d „ 3AM d _ 

T =- Q ©-z 23 

olr 2Mdr 4 dr v ; 

in the phase of broken symmetry, and 

dz A d 



T 



dr AM dr 



CD -z (24) 



in the symmetric phase. The right-hand sides of these equations turn out to depend only 
on z and A, hence the change of z with r (and consequently the change of A with r) is 
uniquely determined for every point on the characteristic curve in the z-A-plane. As for 
the direction of the flow, the right-hand side of Eq. (23) is manifestly negative, hence in 
the broken-symmetry phase the system flows from z — > oo (at r = 0) towards smaller z 
until it arrives, at r = T c , at z — 0. 

The right-hand side of Eq. (24), however, can have both signs, depending on the value 
of z (A is a well-defined function of z for a fixed characteristic curve). It is positive for 
small z and negative for large z, hence the flow with growing temperature is towards an 
intermediate (quasi-) fixed point determined by 

which defines the asymptotic high-temperature values of z and A. For systems which 
enter the symmetric phase after a phase transition at finite temperature, the flow in the 
symmetric phase starts at z = (corresponding to r = T c ) and tends towards larger 
values of z. For systems which are in the symmetric phase at r = 0, the flow starts at 
z — > oo and is directed towards smaller values of z. In both cases, for r — > oo, the fixed 
point (25) is reached asymptotically. 

The characteristic curves for the initial values A = 0.05, 0.5, 2.0 are plotted in Fig. 6, 
together with the line of (quasi-) fixed points (dashed). The line for the validity of the 
HTL approximation is also represented in Fig. 6 (dotted). It intersects every characteristic 
curve once, and at the intersection Eq. (21) is fulfilled. As for the curves corresponding 
to dimensional reduction in Figs. 3-5, they simply correspond to vertical lines at z = z 
(zq = 0.3 throughout most of the paper). The HTL approximation is adequate for the 
part of the characteristic curve to the right of the HTL line (far from the phase transition 
where the fluctuations in A are important), while dimensional reduction applies to the left 
of z — 0.3 (close to the phase transition). 
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0.15 
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Figure 6: The characteristic curves in the z-A-plane for initial values Ao = 0.05, 0.5, 2.0 
(solid curves), together with the line of quasi-fixed points (dashed line) and the HTL 
curve (dotted line) discussed in the text. 



For small z, all loop integrals can be evaluated analytically [7] . One then has that all 
characteristic curves are given by 

X(z) = ^z + 0(z 2 \nz), (26) 

independently of Ao- The universal slope is a reflection of universality at the critical point 
and originates exclusively from the thermal fluctuations. The HTL curve is given by 

\{z) = — + 0{z 2 In z) , (27) 




where we have considered a general v instead of v§ = 0.1 in Eq. (20). It is then clear that 
for sufficiently small z, all characteristic curves lie above the HTL line (for any choice of is), 
so that the HTL approximation does not apply. Finally, for the line of high-temperature 
fixed points, we have 

X(z) = 24z 2 + 0{z 3 ) , (28) 

so that the fixed point line lies below the HTL curve (again for any choice of v) for 
sufficiently small X(z). As a consequence, at high temperatures the HTL approximation 
is satisfactory as long as Ao (and hence A(r) for r — > oo) is sufficiently small. 

The representation of the flow in the z-A-plane makes the high-temperature limit 
particularly transparent through the line of (quasi-) fixed points. First of all, it is clear 
that there are really asymptotic values for A(r) and for M(r)/r. Furthermore, these 
values only depend on Ao and not on the zero temperature mass M , nor on the phase 
the system is in at r = 0. The independence of the physics at high temperatures on 
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M is part of the usual folklore. However, it is also apparent in this representation 
that dimensional reduction and the HTL approximation do not automatically apply at 
high enough temperatures, but that they require the smallness of the coupling constant. 
Finally, the asymptotic value A(r) for r — > oo is always smaller than Ao and can be 
calculated for a given Ao using Eq. (25) and the integrated flow equation for A. In an 
expansion in A around A = 0, one has 



Aoo = lim A(t) = A 




+ 0(A 2 lnA o ) 



(29) 



Analogously, one can determine the asymptotic value for z(t), t — > oo, and obtains for 
the mass in a small-Ao expansion 



lim M 2 (r) = ^l! 




+ o(\ 2 j 



Aor^ 
24 




0(A 2 lnA o ). (30) 



In the HTL approach, Eq. (29) is equally valid, which is a direct consequence of the 
fact that Eqs. (18) and (19) coincide to the order considered here. However, and somewhat 
surprisingly, Parwani [2] finds for the thermal mass 



lim Mi 

T^OO 



iitl( t ) — 



Apr 2 
24 




+ 0(A 2 lnA o ), 



(31) 



which is clearly different from our result (30). In order to decide which of the two re- 
summation schemes is more appropriate in the sense of a faster convergence towards the 
exact result, one would have to calculate the higher-order contributions in both schemes. 
We will not do that here, but rather consider the analogous problem in the limit of large 
N in the O(N) model, which will shed at least some light on this issue. 

In the O(N) model in the symmetric phase, the flow equations take a form nearly 
identical with Eqs. (10, 11), one only has to replace on the right-hand sides 



© 



iV + 2 



N + 8 
9 



(32) 



Repeating the steps leading to Eq. (30) gives the following result for the expansion in Ao 
of the thermal mass at high temperatures, 



lim M 2 (t) 

T— >CX> 



N + 2 A r 2 



24 



1 - 



N 



5\ 3 N + 2 A 
N + 2J if 3 24 



+ 0(A 2 lnA o ). 



(33) 



For iV = 1, we of course recover the result (30). Now, performing the analogue of the 
HTL analysis for the O(N) model, one arrives at 



lim M£ TL (r) 



N + 2 A r 2 



24 



_ i N + 2 *° 

ttV 3 24 



+ O(A 2 lnA ). 



(34) 
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These new results only coincide in the limit iV — * oo. This is a particularly interesting 
limit insofar as the theory is exactly solvable for iV — > oo. The limit has to be taken in 
such a way that 

Ao = Ao (35) 

tend towards a (non-vanishing) constant. As we shall soon see, the expression found 
above for lim r ^ 00 M 2 (r) in the A — * oo limit is correct (to this order in A ). 

As we mentioned above, the 0(A) model can be exactly solved in the large- A limit. 
The result for M(t) and A(r) defined as in Eqs. (3, 4) and (35), is in the symmetric phase 

M 2 (r) = M^ + ^-®'(M(r),r), (36) 

A(r) = -j j-^ , (37) 

I - + -©'(M(r),r) 

where we still have neglected any logarithmic term. We can describe a phase transition at 
a finite temperature by choosing the parameter Mq negative (but the above equations are 
only valid above the critical temperature). Eq. (36) is an implicit equation for M(r). An 
iterative solution leads to the set of the so-called superdaisy diagrams [9]. Eqs. (33) (for 
A — > oo) and (34) are the result of the second iteration in the limit of high temperatures 
and correspond to the daisy diagrams. Eq. (37), when expanded in a geometric series, 
generates the complete series of bubble or chain diagrams where every bubble in itself 
corresponds to the sum of all superdaisy diagrams (with two external vertices). 

It is then easy to verify that the total r-derivative of Eqs. (36) and (37) leads exactly 
to the flow equations (10, 11) with the substitutions (32) in the limit of large A. In other 
words, the integration of the one-loop flow equations leads to the exact solution in the 
large- A limit. This fact gives a lot of confidence in the results of the integration of the 
flow equations also for finite A, in particular, in Eq. (30) for N = 1. Also observe that in 
the HTL formalism, the complete (resummed) perturbative series has to be summed up 
to obtain the results (36) and (37). As a last remark, the use of the total instead of the 
partial r-derivative in the flow equation (10) for the mass, which might have been naively 
expected, can now clearly be seen to be inconsistent, since such flow equations would 
lead to a result different from (36) and (37) in the large- A limit, through overcounting 
of diagrams. 

In large parts of the foregoing discussion, we have neglected the logarithmic contribu- 
tions with the argument that they are negligible for all but extremely high temperatures. 
This approximation, apart from bringing out properties of the flow equations that would 
otherwise hardly be visible in an analytic approach, also simplifies the comparison with 
the HTL approximation and the superdaisy diagram summation a la Dolan-Jackiw [9]. 
For the rest of this section, we will discuss the effects of including the logarithms. 

Beginning with the large- A limit, we reintroduce the logarithms in Eqs. (36) and 
(37) simply by omitting the primes on the loop integrals. These integrals are then diver- 
gent, and we renormalize by subtracting the pole terms as in a naive MS renormalization 
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Figure 7: Solutions of Eqs. (38, 39) for the temperature dependence of the coupling 
constant A and z(t) = M(r)/r in the limit of large N. At temperature r = 0, the system 
is in the symmetric phase with the initial conditions Ao = 2.0 and M = fi. From Eqs. 
(38, 39) one then finds that A r = A = 2.0 and M 2 = [1 + l/(4vr) 2 ] . However, for the 
plot we have neglected M 2 in Eq. (38) since it does not contribute in the limit of high 
temperature which is our primary concern here. 



scheme with temperature-dependent counterterms (we may also subtract the correspond- 
ing ln(47r) — 7^, as in the MS scheme). The result is 

(D'(M(r),r) + Ar2[T 



M 2 (r) 
A(r) 



(4tt) ; 



'in _ i 



1 1 

K + 2 



D'(M(r),r)- 



1 



(47T) 



(38) 
(39) 



where M 2 and A r are the renormalized parameters determined by the choice of initial 
conditions Mq and A , and /i sets a (arbitrary) scale for M(r). Taking again the total r- 
derivatives of these equations, one recovers the flow equations as before, but now including 
all logarithmic terms. 

In practice, the difference between Eqs. (38, 39) and Eqs. (36, 37) is only visible for 
very large temperatures. We present the solution of Eqs. (38, 39) for the initial condition 
A = 2.0 in Fig. 7. Notably, the values of A and z do not vary appreciably over a large 
range of temperatures and stay close to the quasi-fixed point determined from Eqs. (36, 
37), i.e., by neglecting the logarithmic terms (for A = 2.0, the quasi-fixed point is at 
z = 0.2570). This is the significance of the quasi-fixed point for the complete system of 
equations (with logarithmic terms). Finally, at ln(r//x) = 78.47, A diverges while z tends 
towards the finite value 2.05. This behavior is a reflection of the Landau pole at high 
momentum in the zero-temperature theory. 

We mention that the logarithmic terms originating from the correct renormalization 
of the corresponding diagrams, also appear in the original work of Dolan and Jackiw [9] 
and in the HTL expansion [2] . However, since the latter expansion is organized in powers 
of v^Ao and its logarithm, the logarithms of the one-loop diagrams are grouped together 
with contributions from two-loop diagrams. 
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Figure 8: Solutions of Eqs. (10, 11) for the temperature dependence of the coupling 
constant A and z(t) = M(r)/r for N = 1. At temperature r = 0, the system is in the 
symmetric phase with the initial conditions Aq = 2.0 and M = 1 (in arbitrary units). 



Turning now to N — 1, we show the solutions of Eqs. (10, 11) over a large range of 
temperatures in Fig. 8, again for Ao = 2.0. Observe that the qualitative behaviour is as 
in Fig. 7, i.e., for large N (note that for Fig. 7 we have neglected in Eq. (38) which is 
why the curves for z(t) look different for small r). There is a plateau of nearly constant 
values of A and z for a large temperature range, close to the quasi-fixed point of Eq. (25) 
(neglecting the logarithmic contributions). For the present case with Ao = 2.0, the quasi- 
fixed point is at z = 0.232. For smaller initial couplings, the plateau is considerably more 
pronounced — flatter as well as wider, and closer to the quasi-fixed point values. For 
very large temperatures, eventually A runs into the analogue of a Landau pole, at a finite 
value of z. We remark that the integration of the differential equations (10, 11) over such 
a large range of temperatures is a highly non-trivial task from a numerical point of view, 
and our confidence in the results obtained is partly based on the very similar behaviour 
of the solutions of the algebraic equations (38, 39), which we checked for different values 
of the initial coupling A . 

On the analytical side, it is easy to see that characteristic curves in the z-A-plane 
exist even after including the logarithmic term. They are determined by the differential 
equation 



dX 3^ 2 d ^ 

dz 2 dz 



(47T) 



:A 2 



1 1 

z dz 



(40) 



where r(dz/dr) has to be taken from Eq. (23) or (24), according to the phase. However, 
analytic solutions of this equation are not feasible. Qualitatively, the solutions approach 
the line of former high-temperature fixed points from both sides at temperatures of the 
order of several M /\/X~o, and move away from the corresponding quasi-fixed point only 
at very much higher temperatures, following the line of quasi-fixed points very closely 
and slowly towards larger values of z and A. They eventually lift off the fixed point line 
to approach the "Landau pole" in A at finite z. 
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The mere existence of the characteristic curves and the quasi-fixed points implies that 
the physics at sufficiently high temperatures becomes, to a very good approximation, 
independent of the initial value of the mass M and the phase at zero temperature, and 
only depends on the value of Ao- The HTL approximation eventually breaks down for 
any value of A because it is not capable of reproducing a Landau pole without further 
resummation. 

4 Conclusions 

In this paper we have considered the different regimes in which a dimensional reduction 
may occur at finite temperature and compared the validity of different computational 
schemes for accessing these domains. We emphasized that there are two distinct regimes 
in which dimensional reduction may occur at high temperature — with the temperature- 
dependent mass much smaller than and of the same order as the characteristic scale gT, 
respectively — pointing out that the latter is non-universal (NHT regime) whereas the 
former is universal (UHT regime), being independent of the zero temperature parameters. 
NHT dimensional reduction can be successfully described using the HTL approximation as 
in this case the leading order infrared divergences can be resummed leading to a thermal 
mass that acts as an infrared cutoff for fluctuations. The resulting behaviour is then 
mean-field like in that no universal fluctuations lead to strong corrections, as in critical 
phenomena. On the contrary, in the UHT regime HTL-type approximations are invalid 
and other methodologies must be sought. We have shown that environmentally friendly 
RGs allow for a complete consistent analysis of both the NHT and UHT regimes. 

The region in which the HTL approximation breaks down is a window around the 
critical temperature whose size monotonically increases as a function of the zero temper- 
ature coupling. Interestingly, we found that for sufficiently large coupling constants, this 
"window" extends to arbitrarily high temperatures, so that, contrary to common folklore, 
the HTL approximation is not guaranteed to be applicable to this regime. The reason for 
the breakdown of the HTL approximation are strong thermal corrections to the coupling 
constant which cannot be taken into account correctly without a further resummation. 
At extremely high temperatures, we furthermore found a Landau pole for the coupling 
constant, which cannot be described by the HTL expansion to any finite order. 

A similar window around the critical temperature indicates the realization of dimen- 
sional reduction. While for small couplings dimensional reduction applies for arbitrarily 
high temperatures, for sufficiently large couplings it is only valid over a finite tempera- 
ture range. This failure of dimensional reduction for high temperatures which seems to 
go against intuition, has its origin in the linear rise of the thermal mass with tempera- 
ture (over a wide range of high temperatures) with a proportionality constant that grows 
monotonically with the zero-temperature coupling. 

One of the most interesting conclusions we can glean from our results is that whether 
or not there is a regime of non-universal dimensional reduction which can be described 
through a HTL-type analysis depends sensitively on the magnitude of the zero-temper- 
ature coupling constant. We saw that for sufficiently strong couplings the HTL window 
around the critical temperature described above, where the HTL approximation is not 
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valid, completely contains the dimensional reduction window where this reduction obtains. 
Just how strong these couplings have to be was seen to depend on how rigorous our 
criteria were for dimensional reduction and validity of the HTL approximation. This is 
an important finding in view of methods that work with dimensionally reduced theories 
and therefore depend on the existence of a regime where the four-dimensional parameters 
can be mapped onto three-dimensional ones through an analytical, HTL-like procedure. 
In the particular case of the electroweak phase transition, the existence of such a regime in 
the broken-symmetry phase seems to be guaranteed [4]. At any rate, our results further 
confirm that an environmentally friendly RG offers a general methodology for considering 
high temperature field theory that does not suffer from the defects of other methodologies. 

One other defect of the HTL approximation vis a vis environmentally friendly RGs is 
in the case where there are Goldstone bosons in the broken phase as these lead to universal 
infrared divergences that cannot be resummed a la HTL approximation. Environmentally 
friendly renormalization in distinction is capable of treating the Goldstone bosons, their 
effect being to lead to a fixed point for the coupling constant other than the Wilson-Fisher 
fixed point [10]. 
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